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Abstract 


Although the construction of composition profiles from x-ray intensity 
bands is an important and reliable technique, it has not been routinely 
practical because it required hours of computer time and complicated 
optimization routines. The Intensity-Band-to-Composition-Prof ile- 
Trans format Ion proposed here utilizes a solution that can be evaluated even 
with a hand-held calculator. This technique Is applicable to thin films 
and thick specimens for which the variation of lattice parameters, linear 
absorption coefficient and reflectivity with composition are known. A 
deconvolution scheme including corrections for the instrumental broadening 
and a Kordoublet is discussed. 
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Introduction 


The x-ray diffraction intensity band is a compos it ionally broadened 
diffraction peak. The ideal diffraction peak is sharp, narrow, and sym- 
metric. Such a peak is observed by diffraction from a homogeneous phase 
in which the inter-atomic spaclngs are uniform. The intensity band is 
observed by diffraction from a nonhomogeneous phase in which the variations 
in composition result in a range of inter-atomic spaclngs. 

The basic relationships between the x-ray diffraction band and the 
composition-depth profile have been understood for some time. An intrigu- 
ing aspect of these relationships is the hyper- sensitivity of the intensity 
band to the shape of the profile. A number of techniques have been 
developed by investigators seeking to use this sensitivity to construct 
high precision profiles (1-3). The difficulties encountered have been 
twofold: complications due to intensity broadening, and prohibitive 

computational requirements* The seriousness of these difficulties has 
restricted the use of these techniques to academic pursuits. 

The entire breadth of the intensity band is not due solely to composi- 
tional broadening. Instrumental broadening and specimen broadening also 
contribute to the intensity band. To Interpret the Intensity band cor- 
rectly this non-compo8itional broadening must be accounted for. In the 
past simulation techniques have most successfully handled this problem. 

The compositional broadening is calculated from an assumed composition- 
depth profile and convoluted with an analytic broadening function. The 
assumed composition-depth profile is then varied by trial and error until 
there is adequate agreement between the simulation and the experimental 
Intensities (2). 

The simulation techniques have also been the most accurate interpre- 
ters of the intensity band. However, the complexity of the calculations 
and its trial and error dependence necessitated lengthy calculations that 
even with high-speed computers and non-linear regression have been prohibi- 
tively long (3). 

The technique presented here has simple computational requirements and 
is as accurate and flexible as the simulation techniques. It first removes 
the specimen and instrumental broadening by deconvolution and then trans- 
forms the deconvoluted intensities into a composition-depth profile using 
the relationships developed In the theory section of this paper. 


Deconvolution 


Theory 


In a diffractometer the diffracted x-ray Intensity is measured as the 
number of x-ray quanta per second incident on the detector, or, when nearly 
monochromatic x-rays are being used, It is more conveniently expressed as 
power. This power has been broadened by various sources from the instru- 
ment and sample. We are interested In Isolating the compositional broaden- 
ing from all other sources. 

The deconvolution technique Is predicated on expressing the broadened 
power as a summation. To do so we must first adopt two artificial con- 
cepts. The first is that the broadened power, which has been collected at 
arbitrary 28 positions, represents the power-per-degree 26 as if it had 
been collected with an Infinitely narrow receiver slit. The second concept 
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is that the unbroadened power has the form of an array of contiguous 
channels of arbitrary widths and that the power-per-degree within each 
channel is constant. 

Using these concepts we can consider the following convolution 
equation: 



dP b 

where d20 18 the broadened power-per-degree 26, g is a normalized dis- 

ijp 

tributlon function, ^e un broadened power-per-degree 20, and e 

is an error term that is associated primarily with counting statistics. 
Because the form of the unbroadened power is an array of contiguous 
channels, the broadened power may also be expressed by: 



dP 

Where d20 18 the P ower “P eiw degree 20 for channel m and G nm is a 

dimensionless quantity given by: 
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r b i 
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•^—00 1 


J rb 

I “ g(20 - t)dt 
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where the limits of integration, b m , are the channel boundaries. 
It can be shown from Eq. (2) that for the Lagrange function 


F "S i + £ y s 


n-1 


s-1 


(3) 


( 4 ) 


where y 8 is a Lagrange multiplier and $ s is an arbitrary constraint, 
reasonable values for are given by: 
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( 5 ) 


( dP \* 

is a suitable estimate of the unbroadened power, F* ^ are 

cne rirst partlals of the Lagrange function with respect to the unbroadened 
power-per-degree 20, and [F"]^ is the natrlx inverse of the array of 
second partlals of the Lagrange function with respect to unbroadened power- 
per-degree 20. 

When applied In an Iterative fashion Eq. (5) will eventually yield 
reasonable values for the unbroadened power provided that the Lagrange 
multipliers, u s , and the arbitrary constraints, <j» a , are properly chosen. 


It is necessary that at least one constraint be employed to keep the 
N 2 

sum of errors squared, £ e*, from vanishing. If this happens the 
n-1 

unbroadened power-per-degree 20 becomes dominated by large spurious 
oscillations. We suggest that an appropriate constraint would be: 


* ■ (s (Hf)' * ! ■ ° 0 
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where R is the square of the n-fold distance between the successive 
iterations of Eq. (5). This constraint causes the convergence of Eq. (5) 
to progress uniformly. Suitable Lagrange multipliers for this constraint 
lie between 0.1 and 1. 


If the normalized distribution g is the non-compositional broadening 
measured at the intensity-band extremes or from homogeneous phases, then 
the unbroadened power given by Eq. (5) has only compositional broadening 
remaining. The normalized distribution g is not restricted to any 
particular analytic form and, hence, can be chosen to fit any line-shape, 
singlet or doublet. 


Intensity Band to Composition-Depth Profile Transformation 

The technique employed here is based on the well-known equation for 
the power diffracted by an infinitesimally thin volume element at a 
depth x: 


dP - 2P o «X -1 exp(r^)Q dx 


(7) 


where dP is the diffracted power (energy 
incident beam power, 6 is the d-spacing, 
term, X is the x-ray wavelength, n x ia 
the thickness of the volume element. 


per unit time) , P is the 
Q is a pseudo-reflectivity 
an absorption term, and dx 


is 
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The pseudo-reflectivity term, Q, for a polycrystalline material and a 
monochromator with a Bragg angle of 2a is given by: 


Q 


r 2 2 

e 1 4- cos 2a cos 29 

V 2 Sin 26(1 + cos 2 2o) 


F p exp(-2M) 


( 8 ) 


2 *"26 2 

where r fi - 7.94 x 10 cm , V is the volume of the unit cell, F is 
the structure factor, p is the multiplicity, and the exponential term is 
the Debye-Waller thermal factor. 


The absorption term is given by: 





p dx + p*t* 


) 


(9) 


where p is the linear absorption coefficient and t* is the thickness of 
the over layer, if present, and p* is the linear absorption coefficient 
for the overlayer. The first derivative of n with respect to 6 is 
given by: * 

ir" _4ll5x " 1 (1 + v(af) < 10 > 

where the d-spacing gradient term y x is given by: 



Multiplying by the integration factor [2pQ _1 (l + y )] , Eq. (7) 
becomes: x 


( X + Y x ) <35 - -P o expCn^) dn (12) 

We shall define a modified diffracted power l x as the integration of 
either side of Eq. (12) over the range between $ 0 , the d-spacing of the 
surface and 6 X , the d-spacing of the material at the depth x. Thus: 

f 6 x 

“ [ ( S ) 2 *" 1(1 + di < 13 > 

o 

and 

“ P o [ ex P(n o ) - exp(n x )] (14) 
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where 


-1 * * 

n o ■ " 46 o X u t 


( 15 ) 


The quantity J is equal to the counts (x-ray quanta) per second 

measured by the detector. All of the quantities in Eqs. (13) and (14) are 
explicitly known except for x and y x . however , is very small when 

t* - 0 and can be approximated by zero to a very high precision. Thus, 
solving for x in Eq. (14): 


* 1 Ij ["AW)” 1 ta(l - expC-n^P^ 1 I^)] d« (16) 


Usually it is not convenient to determine P Q experimentally. A 
numerical approximation consistent with Eq. (14) can be obtained from 



P o “ exp(-n 0 ) (17) 

because exp( q^) - 0. is the integration of Eq. (13) over the entire 

intensity band— i.e. , it is the total modified power. 

For thin film P Q must be calculated in an iterative fashion from the 
expression: 


p 0 " *t t ex P<no> _ex P(n t )] -1 (18) 

where is the total modified power analogous to B w 

n t " (<!»> t * P*t*) (19) 


t is the film thickness, <p> is the average linear absorption coeffi- 
cient, and $ t is the d-spacing corresponding to the backside of the film. 
Initially a value is assumed for <y>, P Q is calculated from Eq. (18), 
and the profile evaluated using Eq. (16). Subsequently <y> is calculated 
by: 


<y> - T f y dx (20) 


and P Q and the profile are reevaluated. The average linear absorption 
coefficient, P Q , and the profile are evaluated iteratively until the 
change in the profile between successive iterations is satisfactorily 
small (5). 
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Discussion 


Consider a flat sample containing a nonhomogene ous phase such as the 
one depicted In Figure 1. In this example the nonhomgeneous phase Is the 
substrate and a second phase Is an overlayer. The composition gradient Is 
normal to the free-surface and the substrate Is, for x-ray purposes. 
Infinitely thick. 




Depth 


a. Overlayer 

b. Near surface ) Substrate - single phase with 

c. Shallow interior ( a composition gradient 

d. Deep interior j 

Figure 1. A specimen with a composition gradient 
in the substrate. The overlayer and substrate are 
different phases. 


When the incident beam enters the sample, as depicted in Figure 2, no 
diffraction pertinent to the intensity band occurs in the overlayer. The 
sole effect of the overlayer Is absorption. Diffraction occurs in the 
substrate. X-rays diffracted near the surface undergo very little absorp- 
tion and therefore the diffracted intensity from the near surface is 
strong. Diffraction from the interior is weaker because both the incident 
beam and the diffracted beams are heavily absorbed - diffraction from the 
deep interior is completely lost. The angle of diffraction is determined 
by the composition. Each point along the composition gradient will 
diffract the x-rays at a slightly different angle; the combined effect 
is the intensity band. 



Figure 2. 
gradient. 


a. Overlayer 

b. Near surface 

c. Shallow interior; 

d. Deep interior 


no diffraction 
strong diffraction 
weak diffraction 
completely absorbed diffraction 


Diffraction from a specimen with a composition 
Change of 29 with depth is continuous. 
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If we assume that the composition gradient Is normal to the specimen 
surface and that the composition profile Is monotonlc then the transforma- 
tion outlined in the theory section of this paper will construct the com- 
position-depth profile from the intensity band. If the phase being 
examined is in the form of a thin film it must be uniformly thick. The 
accuracy of the transformation will then depend on how well the thickness 
is known. A phase is considered to be in the form of a thin film for jc-ray 
purposes if there is perceptible diffracted intensity from all parts of the 
phase. 

Because this is not an elemental technique the results will be 
ambiguous except for binary composition depth profiles. Quasi-binary 
scenarios such as the dissociation and absorption of a gas by a metal alloy 
are also unambiguous, but cases where there are several diffusing species 
can only be analyzed under special circumstances which are beyond the scope 
of this paper. 

Currently this technique is being employed in diffusion, kinetic 
modeling, and phase transformation studies. Its potential use is much 
broader including semiconductor and quality control applications because of 
its high resolution and its ability to discriminate phases. 

The working depth can be varied by changing the wavelength of the 
incident beam. For most materials the working depth can be varied from a 
few hundred angstroms to several millimeters. This is determined by the 
linear absorption coefficient and is approximately equal to (2y)~ . 

Duplicate specimens of oxidized titanium were analyzed and the results 
of the analysis were compared. Figure 3 presents the composition-depth 
profiles of three sets of duplicate specimens. Specimens in a duplicate 
set were prepared independently. Differences between composition-depth 
profiles for duplicate specimens reflect experimental error as well as the 
inherent precision of the technique. 



transformation. Use of replicates is, for the first time, 
practical because of the ease of data acquisition and analysis* 
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In order to describe the effect of counting statistics on the composi- 
tion depth profile consider again the expression for the modified dif- 
fracted power, Eq. (13). Because the standard deviation of L counts Is 
counts (assuming a Poisson distribution), 



(L ± 1.65 ^)/ T 


2yQ _1 (l + Y x > d« 


( 21 ) 


where 1^ are the 90% confidence limits for I x , L is the number of 
counts, and y is the counting time. The 90% confidence limits on the 
depth are given by: 

X* - jf w" 1 [l- exp(-n 0 ) p” 1 i*JJd« (22) 

where x* are the 90% confidence limits on the depth. An example calcu- 
lation of x* is presented in Figure 4. Errors in the profile due to 
counting statistics tend to be small, being the smallest at x « 0 and 
increasing with depth. 



Figure 4. Confidence limits calculation for an intensity band 
transformation. It is possible to achieve good precision in the 
composition-depth profile using continuous scan data — opposed 
to point counting — as demonstrated by the transformation of an 
intensity band collected at a scan speed of 0.1° per minute 
over 1.5°. 
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The diffractometer used should have moderate to high resolution. 
Diffracted beam monochromators that can separate the Ka^ component are 
desirable. It is also good to choose the narrowest receiver slit available 
and in general to do all things that aid resolution. 

The use of high angle peaks and high resolution diffractometers will 
help preserve the detail contained in the intensity band. The more non- 
compositional broadening present, the more one will need to rely upon the 
deconvolution routine to construct meaningful profiles. Failure to decon- 
volute adequately will produce errors at the extreme compositions as well 
as masking the detailed structure of the composition-depth profile. 

Figure 5 is a composite plot of the transformation of an intensity band 
before and after deconvolution. The transformation of the Intensity band 
before deconvolution has errors at the extreme compositions. 



By making a series of approximations and simplifications, including 
Q - constant and y » constant, it is possible to write: 


_-sin 6n _ 

■ —s 5 1 xn 

2y 


‘- e (& We (&)“ 3 

i -1 ' / j -1 J 


(23) 


where ^ is the depth 


•(4 


is the counts per second at the detector 


arranged in the order such that 


(4 


corresponds to 6 




responds to This expression is similar in several points to that of 

Pines ( 6 ). 


This expression is sufficiently simple so that it could be evaluated 
for small data sets with a hand-held calculator. Assuming that non- 
compositional broadening is negligible we may skip the deconvolution step. 
An example calculation using Eq. (23) is given in Table I and the results 
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are plotted in Figure 6 against the more complete treatment using decon- 
volution and Eq. (16). The agreement is excellent except for a small depth 
(less than 0.3 ym) where the neglecting of instrumental broadening has a 
significant effect. 


Table I. An Example Calculation of Composition Depths Profile 
Using Eq. (23). 


2®i 

/dP\ 

V d5 A 

n 



V im) 

Comp, at Z 

37.47 

0.66 

2.3983 

-0.0022 

-0.0014 

0.006 

25.93 

37.54 


2.3939 

-.0040 

-.0057 

.026 

23.81 

37.60 

2.03 

2.3903 

-.0037 

-.0131 

.061 

21.89 

37.66 

3.88 

2.3866 

-.0040 

-.0285 

.135 

19.90 

37.73 

6.81 

2.3823 

-.0037 

-.0534 

.262 

17.52 

37.78 

9.28 

2.3793 

-.0030 

-.0816 

.417 

15.81 

37.83 

10.01 

2.3763 

-.0039 

-.1209 

.661 

14.12 

37.91 

10.13 

2.3714 

-.0042 

-.1636 

.972 

11.47 

37.97 

9.51 

2.3678 

-.0048 

-.2093 

1.377 

9.56 

38.07 

8.50 

2.3618 


-.2651 

2.046 

6.63 

38.19 

7.55 

2.3547 

HSS 

-.3122 

2.912 

3.65 

38.28 

6.93 

2.3494 


-.3450 

3.922 

1.92 

38.35 

6.19 

2.3452 


-.3687 

5.331 

.93 

38.41 

4.91 

2.3417 

-.0029 

-.3829 

7.574 

.36 

38.44 

2.06 

2.3394 

-.0023 

-.3876 

10.482 

.13 

38.48 

.34 

2.3372 

-.0025 

-.3884 

12.888 

.04 

38.53 

.22 

2.3344 


-.3887 

00 

- 


30 


t 



0 1 2 3 4 5 6 


Depth (microns) 

Figure 6. Comparison of simple power summation transformation with 
the complete technique. The simple power summation transformation 
is useful for checking and quick analysis. The agreement is largely 
fortuitous because the various errors tend to cancel in this example. 
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Concluding Remarks 


Obtaining the composition profiles from x-ray diffraction intensity 
bands is easy, both experimentally and analytically. Though the technique 
is not elemental and hence is usually restricted to binary or quasi-binary 
P r °fil es » the inherent discrimination of phases enables observation of 
clear composition-depth profiles when overlayers or sublayers are present. 
And finally, the composition-depth profiles found by this technique are 
very precise in working depths that range from angstroms to millimeters. 
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